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Abstract Recent measurements of the lY-boson 
charge asymmetry and of the Z’-boson production cross 
sections, performed at the Tevatron collider in Run II 
by the DO and CDF collaborations, are studied using 
the HERAFitter framework to assess their impact oir 
the proton parton distribution functions (PDFs). The 
Tevatron measurements, together with deep-inelastic 
scattering data from HERA, are included in a QCD 
analysis performed at next-to-leading order, and com¬ 
pared to the predictions obtained using other PDF sets 
from different groups. Good agreement between meas¬ 
urements and theoretical predictions is observed. The 
Tevatron data provide significant constraints on the d- 
valence quark distribution. 
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1 Introduction 

Accurate knowledge of the parton distribution func¬ 
tions (PDFs) is essential for predictions at hadron col¬ 
liders. The primary source of information on the pro¬ 
ton PDFs comes from deep-inelastic scattering (DIS). 
Measurements at fixed target experiments and at the 
HERA e^p collider provide constraints on the quark 
and gluon densities, and discrimination of the quark 
flavours. The DIS proton data mostly constrain the u- 
type quark density, due to the greater couplings to the 
photon at low absolute four momentum transfers, Q^, 
whereas the d-type quark densities are only constrained 
at high with limited precision. Even more chal¬ 
lenging is the separation of the d-valence quark dens¬ 
ity which relies on the HERA e~^ charge current data, 
which are statistically limited in the published HERA I 
combined data [1] . A better flavour separation is needed 
to challenge the limits of precision physics at the LHC. 

Drell-Yan production of W and Z bosons in proton- 
antiproton and proton-proton collisions can provide ad¬ 
ditional information on the d-quark PDFs. At lead¬ 
ing order (LO) in QCD, the Drell-Yan processes probe 
the PDFs at energy scales Q corresponding to the bo¬ 
son masses, my = rnw and my = rnz, and mo¬ 
mentum fractions carried by the interacting partons of 
2 ^ 1,2 = my /, where \fS is the centre-of-mass en¬ 
ergy and y is the boson rapidity. 

At the Tevatron proton-antiproton collider, the pro¬ 
duction of W and Z bosons is dominated by valence- 
quark interactions. The Z-boson production has sim¬ 
ilar couplings for uu and dd fusion processes, whereas 
W bosons are produced predominantly by ud and du 
fusions for W~^ and W~ bosons, respectively. Various 
measurements of ^-boson inclusive production and of 
VF-boson charge asymmetry have been reported by the 
DO and CDF collaborations [2-7]. Some of these data 
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samples were included in previous PDF studies [8-11]. 
The addition of the Tevatron data resulted in improved 
PDFs, but some tensions were observed with global 
PDF fits [12-15], 

In this paper the data collected at the Tevatron col¬ 
lider in Run II are analysed to assess their impact on 
the PDFs. The assumptions of the correlation model of 
the experimental systematic uncertainties are revised 
with respect to the recommendation of the experiments, 
leading to improved agreement with the theoretical pre¬ 
dictions. The analysis is performed using the HERA¬ 
Fitter framework [1, 16-18] at next-to-leading order 
(NLO) QCD. The Tevatron W- and Z-boson meas¬ 
urements are also compared to predictions evaluated 
with the recent PDF sets CTlOnlo [8], MMHT2014 [9] 
and NNPDF3.0 [10]. The impact of the Tevatron data 
on PDFs is studied using hessian profiling [19] and 
Bayesian reweighting [20-22] techniques. The profiling 
of PDF uncertainties is generalised to the case of asym¬ 
metric PDF uncertainties. 

This paper is organised as follows: the data samples 
are introduced in Sec. 2 and the theoretical predictions 
are discussed in Sec. 3. The QCD analysis settings and 
the methods for comparing data with predictions based 
on existing PDFs are discussed in Sec. 4. Section 5 re¬ 
ports the results of the PDF analysis. The results ob¬ 
tained in the paper are summarised in Sec. 6. 

2 Experimental Measurements 

2.1 Data Sets 

The most recent measurements of VF-boson charge 
asymmetry and Z-boson inclusive production per¬ 
formed in Run II of the Tevatron collider are considered 
in this study. They include the Z-boson differential 
cross section as a function of rapidity, measured by the 
DO collaboration with 0.4 fb~^ of integrated luminos¬ 
ity in the Z —>■ ee channel [2]; the Z-boson differential 
cross section as a function of rapidity, measured by the 
CDF collaboration with 2.1 fb“^ of integrated lumin¬ 
osity in the Z ^ ee channel [3]; the charge asymmetry 
of muons as a function of rapidity in IT —>■ /ii/ decays, 
measured by the DO collaboration with 7.3 fb“^ of in¬ 
tegrated luminosity [4]; the VF-boson charge asymmetry 
as a function of rapidity in the W ^ ev decay channel, 
measured by the CDF collaboration with 1 fb“^ of in¬ 
tegrated luminosity [5]; the VF-boson charge asymmetry 
as a function of rapidity in the W ^ ev decay chan¬ 
nel, measured by the DO collaboration with 9.7 fb“^ of 
integrated luminosity [6]. These measurements super¬ 
sede the previous Run II Tevatron measurements of VF- 
boson charge asymmetry and Z-boson inclusive produc¬ 


tion. Recently, the DO collaboration has also released a 
measurement of the charge asymmetry of electrons as 
a function of rapidity in IT —>■ ez/ decays [7]. However, 
this measurement is performed with the same data set 
and event selection as the measurement of Ref. [6] , and 
cannot be included simultaneously in a PDF fit without 
provision of the correlation information. The Tevatron 
W- and Z-boson measurements considered in this study 
are summarised in Tab. 1. 

Besides the Tevatron W- and Z-boson measure¬ 
ments, the HERA I combined measurements of the 
inclusive DIS neutral- and charged-current cross sec¬ 
tions measured by the HI and ZEUS experiments [1] 
are used in this study. The neutral-current measure¬ 
ments cover a wide range in Bjorken-cc and Q^, which 
is essential for the determination of PDFs, whereas the 
charged-current measurements provide further inform¬ 
ation to disentangle the contributions in PDFs from 
u-type and d-type quarks and anti-quarks at a: > 0.01. 
The DIS data are required to be in the kinematic re¬ 
gion = 7.5 GeV^, where perturbative QCD 

calculations are reliable. 

2.2 Experimental Uncertainties 

Statistical uncertainties are considered to be uncorrel¬ 
ated between bins, with the exception of the DO meas¬ 
urement of IF-boson charge asymmetry of Ref. [6], for 
which bin-to-bin statistical correlations are provided. 

In general, the correlation model of the experi¬ 
mental uncertainties recommended by the Tevatron ex¬ 
periments is adapted and followed in the QCD analysis, 
with the exception of the experimental systematic un¬ 
certainties related to trigger and lepton identification 
efficiencies. These uncertainties are provided by the DO 
and CDF experiments in the form of total uncertainties 
in each bin of the measurements. However, the trigger 
and lepton identification corrections are estimated from 
data, and they are influenced, among other effects, by 
statistically uncorrelated bin-to-bin fluctuations. Since 
the exact bin-to-bin correlation pattern of these un¬ 
certainties is not provided by the experiments, a con¬ 
servative approach is followed in this study, and the 
uncertainties related to trigger and lepton identifica¬ 
tion efficiencies are treated as uncorrelated bin-to-bin 
for the nominal fit. According to this prescription, the 
following uncertainties are treated as uncorrelated bin- 
to-bin: the central- and forward-electron identification 
efficiencies of Ref. [3] , the trigger isolation efficiency of 
Ref. [4], the trigger and electron identification efficien¬ 
cies of Ref. [5], and the electron identification, charge 
misidentification and positron to electron efficiency cor¬ 
rections of Ref. [7]. 
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All the other experimental systematic uncertain¬ 
ties are considered fully correlated bin-to-bin, with the 
exception of the DO measurement of fD-boson charge 
asymmetry of Ref. [ 6 ], where the total experimental 
systematic uncertainty is treated as bin-to-bin uncor¬ 
related, as recommended by the DO experiment, and 
for the electron charge asymmetry in W ^ ev decays 
of Ref. [7] , where the uncertainty of the unfolding pro¬ 
cedure due to the limited statistics of the Monte Carlo 
(MC) sample is treated as uncorrelated bin-to-bin. The 
dependence of the measured asymmetry on the PDF set 
used to reconstruct the fF-boson rapidity was studied 
in the DO measurement of W -boson charge asymmetry 
of Ref. [ 6 ]. In this paper, the VF-boson charge asym¬ 
metry extracted with the CTEQ 6.6 PDF set is used 
as the central value, and the 22 CTEQ 6.6 positive and 
negative PDF eigenvector variations are considered as 
bin-to-bin correlated systematic uncertainties. 

For the two measurements of ^-boson differential 
cross section as a function of rapidity, statistical uncer¬ 
tainties are scaled to the expected number of events as¬ 
suming they are Poisson distributed. The experimental 
systematic uncertainties are treated as multiplicative, 
and linearly scaled to the expected cross sections, ex¬ 
cept for the background uncertainties which are treated 
as additive, and are not scaled. For the measurements 
of IF-boson and lepton charge asymmetry, all the un¬ 
certainties are treated as additive, and are not scaled. 

The statistical uncertainties of the HERA I data are 
treated as uncorrelated and scaled to the expected num¬ 
ber of events assuming Poisson distribution, whereas 
the experimental systematic uncertainties are fully cor¬ 
related and are scaled linearly to the expected cross 
sections. 

3 Theoretical Predictions 

The theoretical predictions corresponding to the Tevat¬ 
ron measurements of kF-boson charge asymmetry and 
Z-boson inclusive production are included in the fits 
using APPLGRID [23, 24] files. These predictions have 
been evaluated with MCFM [25, 26] at NLO QCD ac¬ 
cording to the phase-space definitions of each measure¬ 
ment, which are as follows: the DO and CDF measure¬ 
ments of the Z-boson differential cross section as a func¬ 
tion of rapidity are defined in the full kinematic range 
of the decay leptons, without any requirements on the 
rapidity and pt of the leptons. In the DO measurement, 
the invariant mass of the dielectron system is defined 
in the range 71 < niee <111 CeV, whereas in the CDF 
measurement, it is defined in the range 66 < mee < 
116 CeV. The charge asymmetry of muons as a func¬ 
tion of rapidity in W ^ pv decays, and the charge 


asymmetry of electrons in IF —> ei^ decays, measured 
by the DO experiment, are defined with pfp > 25 CeV 
and > 25 CeV. The IF-boson charge asymmetry as 
a function of rapidity in the IF —> decay channel, 

measured by the CDF collaboration, is defined in the 
full kinematic range, without any requirements on the 
lepton rapidity and pr- The corresponding DO measure¬ 
ment of the IF-boson charge asymmetry in the IF —>■ 
decay channel is defined in a kinematic region where 
the charged lepton and the neutrino are required to 
have Pt > 25 CeV without further requirements on the 
lepton rapidity. The kinematic requirements of the Tev¬ 
atron IF- and Z-boson measurements are summarised 
in Tab. 1. Notice that the CDF and DO measurements of 
IF-boson charge asymmetry in the W ^ ev decay chan¬ 
nel of Refs. [5] and [ 6 ] are defined in different kinematic 
regions and they should not be compared without ex¬ 
trapolating them to a common phase space. Tables of 
the Tevatron measurements, with updated correlation 
model, and corresponding APPLCRID theoretical pre¬ 
dictions are publicly available at herafitter.org. 

The QCD predictions for the DIS cross sections 
are evaluated by solving the DCLAP evolution equa¬ 
tions [27-32] at NLO in the MS scheme [33] using 
the QCDNUM program [34] with the renormalisation 
and factorisation scales set to Q^. The light quark 
coefficient functions are calculated in QCDNUM. The 
heavy c- and 5-quark distributions are dynamically 
generated, and the corresponding coefficient functions 
for the neutral-current processes with 7 * exchange are 
calculated in the general-mass variable-flavour-number 
scheme [35-37], with up to five active quark flavours. 
For the charged-current processes and the neutral- 
current processes with a Z contribution, the heavy 
quarks are treated as massless. 


4 QCD Analysis 

4.1 PDF Fit Settings 

The QCD analysis and PDF extraction is performed 
with the open-source framework HERAFitter. The 
charm mass is set to rric = 1.38 CeV, as estimated 
from HERA charm production cross section [38], and 
the bottom mass to mi, = 4.75 CeV [39]. The strong- 
interaction coupling constant at the Z boson mass, 
Ois{Mz), is set to 0.118, and two-loop order is used for 
the running of Ug. 

The PDFs for the gluon, m- valence, d-valence, u, 
d quark densities are parametrised at the input scale 
Qq = 1.7 CeV^ as follows: 
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xf{x) = AfX^^{l — x)^^ {1 + DfX + EfX^)e^^^ ; (1) 

f Uy , dy ^ g ^ d. 

The contribution of the s-quark density is taken to be 
proportional to the d-quark density by setting xs{x) = 
rsxd{x), with rg = 1.0, as suggested in Ref. [40]. The 
strange and anti-strange quark densities are taken to be 
equal: xs{x) = xs(x). The normalisation of the xuy(x) 
(xdy(x)) valence-quark density, Ay^ is determ¬ 

ined by the quark-counting sum rule, whereas the nor¬ 
malisation of the gluon density, Ag, is determined by 
the momentum sum rule. The a; —>■ 0 limit of the u- 
to d- sea quark densities is fixed to unity by setting 
Ey — Rj and Ay — A^. 

A function used for the data to theory compar¬ 
ison is defined as in Ref. [1], with an additional pen¬ 
alty term as described in Ref. [41], and minimised with 
MINUIT [42] to extract the PDFs from the data. 


4.2 PDF Profiling 


where and 12 *^ = ^ 

Fik ) are determined from the shifts of predictions cor¬ 
responding to up and down {F^k~) PDF uncer¬ 

tainty eigenvectors. 

The minimisation of Eq. 2 in its original form leads 
to a system of linear equations. The generalised func¬ 
tion, with asymmetric PDF uncertainties, is minimised 
iteratively: the values of F-^~^ are updated using /3fc,th 
from the previous iteration and following the substitu¬ 
tion of Eq. 3. Several iterations are required to con¬ 
verge, and the procedure is verified using the MINUIT 
program which yields identical results. 

The value at the minimum of the function 
provides a compatibility test of the data and theory. 
In addition, the values at the minimum of the nuis¬ 
ance parameters Fe interpreted as optimisa¬ 

tion (“profiling”) of PDFs to describe the data [19]. 
Explicitly, the profiled central PDF set /q is given by 


/o = /o+E'3“tS 

k 



/□min 

Pfc,th 




(4) 


The impact of a new data set on a given PDF set can be 
quantitatively estimated with a profiling procedure [19]. 
The profiling is performed using a function which 
includes both the experimental uncertainties and the 
theoretical uncertainties arising from PDF variations: 


X (/3expj/3th) — 


^data 

E 




_th 




+ E '^l.exp + E 

3 k 



( 2 ) 


where /o is the original central PDF set and repres¬ 
ents the eigenvector sets corresponding to up and down 
variations. 

The shifted PDFs have reduced uncertainties. In 
general, the shifted eigenvectors are no longer ortho¬ 
gonal, but can be transformed to an orthogonal repres¬ 
entation using a standard diagonalisation procedure, as 
in Ref. [43]. In this method the covariance matrix C of 
the PDF nuisance parameters is diagonalised as 


= (G'/3th)"’G'/3th = (/3^h)^/3^h : (5) 


The correlated experimental and theoretical uncertain¬ 
ties are included using the nuisance parameter vectors 
/3exp and ,3th, respectively. Their influence on the data 
and theory predictions is described by the and 
matrices. The index i runs over all Adata data points, 
whereas the index j (k) corresponds to the experi¬ 
mental (theoretical) uncertainty nuisance parameters. 
The measurements and the uncorrelated experimental 
uncertainties are given by and Ai , respectively, 
and the theory predictions are af'. The x^ function 
of Eq. 2 can be generalised to account for asymmetric 
PDF uncertainties: 

(3) 


where G is an orthogonal matrix, D is a positive defin¬ 
ite diagonal matrix, and '/D is a diagonal matrix built 
of \/Dii. The matrices G and D can be constructed 
using the eigenvectors and eigenvalues of the matrix 
G. The transformation G' can be adjusted, using or¬ 
thogonal transformations, to keep the new eigenvector 
basis aligned along the original as much as possible. As 
a result of this adjustment, the transformation matrix 
can take a triangular form with all diagonal elements 
greater than zero. 

The method can be extended to PDF sets with 
asymmetric uncertainties: the transformation matrix is 
determined using symmetrised uncertainties as in Eq. 5, 
and the orthogonal up and down PDF eigenvectors ff' 
and /“ are calculated as 


Fik FI^ -\- , 
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ft' = /o + E,- G', 
/■' = /o - E,- G', 



+ GG 




/ + + /--2/o>j 
2 

/++/r- 2 /o\ 


those which do not describe the new data), and they do 
not contribute to the ensemble any longer. The number 
of effective replicas, Ngg, of a reweighted set is quanti¬ 
fied by the Shannon entropy 


4.3 Bayesian Reweighting 


TVeff = exp 


^ A^rep 

— X! ln(iVrep/tCfc) > . 

k=l J 


(9) 


An alternative approach to assess the impact of new 
data on PDFs is the Bayesian reweighting technique, 
first proposed in Ref. [20] and further developed by the 
NNPDF collaboration [21, 22]. The Bayesian reweight¬ 
ing can be applied to PDF sets provided in the form of 
MC replicas, such as the NNPDF3.0 set [10]. Recently, 
a variant of the method which can be used with PDFs 
provided in the eigenvector representation has been de¬ 
veloped [44] and is also available in HERAFitter. 

The Bayesian reweighting is based on the assump¬ 
tion that an ensemble of MC replicas provides a rep¬ 
resentation of the probability distribution in the space 
of PDFs. For a given PDF set with TVi-ep replicas {fk}, 
with k = 1, 2,..., iVi-ep, the central value for a general 
observable, 0{{fk}), is estimated as the average of the 
predictions obtained from the ensemble: 

Nrep 

(^) = E ■ (6) 

fc=l 

With the inclusion of new data, the probability dis¬ 
tribution associated with the original PDF set is modi¬ 
fied according to Bayes Theorem. For each replica k, a 
weight Wk is obtained from the function according 
to: 


Wk = 


[y^)2(-^data l)g 2 Xfc 


(7) 


where A^data is the number of new data points and Xk is 
the x^ value between data and predictions correspond¬ 
ing to the fc-th PDF replica. 

The prediction for a given observable, after the in¬ 
clusion of the new data, is evaluated as the weighted 
average of predictions obtained from the ensemble: 


A^rep 

^ E ■ (8) 

fc=i 

The reweighting procedure is very fast and results 
in a new, updated, MC PDF set. Some of the replicas 
of the PDF set may have very small weights (typically 


An un-weighting procedure can be performed on the 
MC set such that PDFs with small weights are sup¬ 
pressed and a new set is produced, which has unit 
weight for all PDF replicas in addition to statistically 
reproducing the averages from Eq. 8. 


5 Results 

The QCD fit analysis described in Sec. 4 is performed 
on the Tevatron W- and Z-boson data, together with 
the HERA I data. The fit is used to study the com¬ 
patibility of the data with NLO QCD predictions, and 
to assess the impact of the Tevatron data on PDFs. 
The prohling and the reweighting techniques are used 
to asses the impact of the Tevatron data on various 
PDF sets. 

The optimal parametrisation for the PDF fit is 
found through a parametrisation scan, a procedure first 
introduced in Ref. [1] . The scan is performed by starting 
from a parametrisation with a basic polynomial form, 
where Df, Ej, and the exponential parameters F/ of 
Eq. 1 are set to zero. After application of the quark¬ 
counting and momentum sum rules, and of the a; —>■ 0 
constraints on u and d, the initial PDFs parametrisa¬ 
tion has 10 free parameters. The 15 Df, Ef and Ff 
additional parameters are allowed to vary, one para¬ 
meter at a time, and the parameter which induces the 
largest reduction of Xmin added as a free parameter 
for the next iteration of the scan. The PDF fits which 
lead to solutions with negative high-a; PDFs are dis¬ 
carded. For each PDF, the exponential term and 
the polynomial term {l-\-DfX + Efx‘^) are considered as 
mutually exclusive, that is, when the exponential term 
is preferred, the polynomial term is no longer considered 
in the scan, and vice versa. The procedure is stopped 
when the reduction in the xtin value, is less 

than unity. 

Table 2 shows the results of the parametrisation 
study: the parameters which induce the largest ^xtin 
are, in order, Fd„, Fu„, Dg, Dj, and Du. The optimal 
parametrisation found with this procedure has 15 free 
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parameters, and the PDFs are expressed as: 


xg{x) = AgX^o{\ —Dgx)\ ( 10 ) 

xuy{x) = (1 — ; (11) 

xdy{x) = Ad^x^'^'‘ (1 — a;)*"'^” ; ( 12 ) 

xu{x) = AyX^^ {1 — x)^^ {1 + Dyx); (13) 

xd{x) = AdX^^{l — x)^^{l + D^x). (14) 


The parametrisation of Eqs. (10-14) is used for a 
fit to the HERA 1 data, and for a combined fit to the 
HERA I and Tevatron W- and Z-boson data. Table 3 
shows the Xmin P®i' degrees of freedom (dof) of the two 
fits. The contribution to the total Xmin of data set, 
henceforth referred to as partial x^, is also shown. The 
inclusion of the Tevatron W- and Z-boson data in the 
fit, which corresponds to 93 additional points, results 
in an increase of about 110 in the overall Xmin of fde 
fit, and the partial x^ per number of points of each of 
the Tevatron and HERA 1 data set is close to unity. 

Figures 1 and 2 show the Tevatron Z- and IF- 
boson measurements, respectively, compared to the the¬ 
oretical predictions evaluated with the PDFs extracted 
from the combined fit to the HERA 1 and Tevatron 
data. 

The central value and the uncertainties of the PDFs 
are evaluated with MC replicas [45]: the data points 
are smeared using Gaussian distributions, according to 
their experimental uncertainties, and the PDF fit is re¬ 
peated 1000 times, using different random seeds for the 
smearing. The central PDFs are calculated as the aver¬ 
age of the replicas and the PDF uncertainties are calcu¬ 
lated from their standard deviation. Figure 3 shows the 
comparison of the PDFs extracted with the MC-replica 
method by htting the HERA 1 data, and by fitting the 
HERA I and Tevatron W- and Z-boson data. Figure 4 
shows the comparison of the relative uncertainty of the 
two PDFs. A significant reduction of the PDF uncer¬ 
tainties is observed in the fit which includes the Tev¬ 
atron W- and ^-boson measurements, in particular for 
the valence quarks and d quarks. 

A ht of the HERA I and Tevatron W- and Z-boson 
data with the same settings, but with a correlation 
model in which trigger and identification uncertainties 
are treated as correlated bin-to-bin, yields very similar 
central PDFs and PDF uncertainties [not shown]. The 
X^ of all the data sets are also very similar, except for 
the x^ of the CDF IF-boson asymmetry measurements, 
which is about twice as large. 

The IF-boson charge asymmetries rely on the recon¬ 
struction of the VF-boson rapidity, which is measured 
assuming a fixed VF-boson mass, and inferring the un¬ 
measured longitudinal momentum of the neutrino on 
a statistical basis [5, 6 ]. The reconstruction of the IF- 
boson rapidity introduces a model dependence in the 


measurement. To study the possible bias due to the 
IF-boson rapidity reconstruction, an alternative fit is 
performed in which the VF-boson charge asymmetries 
measured by CDF and DO are excluded, and the latest 
DO measurement of the electron asymmetry is included. 
The Xmin / the partial x^ of the Ht are shown 

in Tab. 5. Also for this fit the partial x^ of each of 
the Tevatron and HERA I data set is close to unity. 
The dy PDF determined from the fit is shown in Fig. 5, 
and compared to the nominal fit. The fit to the lepton 
asymmetry data yields very compatible results, but the 
uncertainties on the dy PDF are up to twice as large. 

The impact of posterior inclusion of the Tevatron 
W- and Z-boson measurements on the PDF uncer¬ 
tainties as estimated by CTlOnlo, MMHT2014, and 
NNPDF3.0 is assessed by profiling and reweighting. For 
consistency with the other PDFs, the uncertainties of 
the CTlOnlo PDFs are scaled to 68 % confidence in¬ 
terval by applying a factor of 1.645. The three PDF 
sets already include the CDF and DO Z-boson differ¬ 
ential cross sections as a function of rapidity, and the 
MMHT2014 fit also includes the DO muon charge asym¬ 
metry in IF —>■ decays and the CDF IF charge asym¬ 
metry in the W ^ ev decay channel. Only the measure¬ 
ments that are not included in each of the PDF sets are 
considered for the corresponding profiling or reweight¬ 
ing. The compatibility of the Tevatron data with the 
CTlOnlo, MMHT2014 and NNPDF3.0 sets is tested by 
evaluating the x^ function of Eq. (2), accounting for 
asymmetric PDF uncertainties according to Eq. (3). To 
perform this calculation for the NNPDF3.0 set, the co- 
variance matrix for the predictions is decomposed using 
the eigenvector representation. Table 4 shows the com¬ 
patibility between the Tevatron measurements and the 
above PDF sets, together with the partial x^ of each 
data set. The partial x^ per number of points of each 
of the Tevatron data set, and the total x^ / dof, are 
close to unity for all the PDFs, when the x^ evalu¬ 
ation includes the PDF uncertainties. The quality of 
the agreement significantly deteriorates if the x^ evalu¬ 
ation neglects the PDF uncertainty. This effect is more 
pronounced for the CTlOnlo and NNPDF3.0 sets which 
include fewer data from the Tevatron. This indicates the 
significant constraining power of the Tevatron data. 

The CTlOnlo and MMHT2014 PDFs are profiled 
according to Eq. (4). The results of the profiling on 
the d-valence PDFs, and on their relative uncertainties, 
are shown in Fig. 6 . The profiling affects the shape of 
the distribution more for the CTlOnlo when compared 
to MMHT2014 set. Significant reduction of the uncer¬ 
tainties is observed for both sets, in particular in the 
low- and medium-a; range. The NNPDF3.0 PDFs are 
reweighted to the Tevatron data. The number of effect- 
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ive replicas remaining after reweighting, is only 1 
and hence the resulting PDFs are not shown. 

The original and profiled d-valence PDFs, and the 
result of the fit to the HERA I and Tevatron W- and Z- 
boson data, are compared in Figs. 7 and 8 , respectively. 
The profiling using Tevatron data improves agreement 
of the d-valence distribution between the MMHT2014 
and CTlOnlo PDF sets. 


6 Summary 

The HERAFitter framework is used to perform a QCD 
analysis of the DIS data from HERA, together with W- 
and Z-boson production measurements performed at 
the Tevatron collider in Run H. The correlation model 
of the systematic uncertainties of the Tevatron data 
is investigated, and a modification is proposed which 
accounts for the statistical nature of some of the sys¬ 
tematic uncertainties. The Tevatron and HERA data 
are well described by a NLO fit with 15 free paramet¬ 
ers, with a new parametrisation of PDFs which adds 
to the basic form a combination of linear and exponen¬ 
tial terms. The impact of the Tevatron W- and Z-boson 
measurements is assessed by comparing PDF uncertain¬ 
ties from a £t to the HERA data alone, and a fit to the 
HERA and Tevatron data. A significant reduction of 
the uncertainties is observed in the latter case, for the 
valence quarks and d quarks in particular. 

The Tevatron measurements are also compared to 
predictions evaluated with modern PDF sets, and the 
impact of the data on the PDFs is assessed using pro¬ 
filing and reweighting techniques. The profiling tech¬ 
niques take into account asymmetric PDF uncertain¬ 
ties. A good agreement between measurements and pre¬ 
dictions is observed, if the PDF uncertainties of the 
predictions are taken into account. After the inclusion 
of the Tevatron data, the PDF uncertainties on the d- 
valence quarks are significantly reduced, especially for 
the PDF sets which include only the Z-boson data from 
the Tevatron, and the agreement between the various 
PDF sets is improved. These findings highlight the im¬ 
portance of the Tevatron W- and Z-boson production 
data to constrain d-quark and valence PDFs, and sug¬ 
gest that the data should be used in the future global 
PDF analyses. All the supporting material to allow fits 
of the Tevatron data, including the updated correlation 
model and the grid files for fast theory calculations, are 
publicly available on the web page of the HERAFitter 
project. 
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Table 1 Summary of the Tevatron W- and Z-boson measurements. For each measurement the observable, the experiment, 
the integrated luminosity, the phase-space definition, the inclusion in the nominal fit, and the corresponding Ref. are shown. 

Observable 

Experiment 

Integrated 

luminosity 

Kinematic 

requirements 

Used in the 
nominal fit 

Ref. 

da{Z)/dy 

DO 

0.4 fb-i 

71 < rriee < 111 GeV 

yes 

[2] 

d(j{Z)fdy 

CDF 

2.1 fb-i 

66 < mee < 116 GeV 

yes 

[3] 

muon charge asymmetry in VF — >■ yv 

DO 

7.3 fb-i 

> 25 GeV, > 25 GeV 

yes 

[4] 

electron charge asymmetry inW ^ ev 

DO 

9.7 fb-i 

Ef. > 25 GeV, p^ > 25 GeV 

no 

[7] 

W charge asymmetry mW ^ ev 

CDF 

1.0 fb-i 

none 

yes 

[5] 

W charge asymmetry in W ^ eu 

DO 

9.7 fb-i 

Ef. > 25 GeV, p^ > 25 GeV 

yes 

[6] 


Table 2 Results of the parametrisation study. For each additional free parameter D, E, and F, of the gluon, u, and 

d PDF, the reduction of of a fit to the Tevatron and HERA I data, , is shown. For each of the fit with n free 

parameters, with n = 10, 11,12,13, 14, the largest shown in bold, and the corresponding parameter is added as a 

free parameter for the n + 1-parameters fit. The fits which leads to negative high-x PDFs are shown in parenthesis, and are 
not considered in the parametrisation study. 


n 

Xmin/dof 

10 

714/633 

11 

654/632 

12 

619/631 

13 

610/630 

14 

607/629 

Free parameter 



■ 

A, min 



dv 

D 

39 

- 

- 

- 

- 

E 

40 

- 

- 

- 

- 

F 

41 

- 

- 

- 

- 

U-u 

D 

< 1 

41 

- 

- 

- 

E 

1 

38 

- 

- 

- 

F 

< 1 

44 

- 

- 

- 

gluon 

D 

25 

16 

9.8 

- 

- 

E 

9 

8 

5.4 

< 0.1 

0.3 

F 

35 

19 

7.2 

- 

- 

d 

D 

2 

< 1 

1.9 

3.0 

- 

E 

(44) 

(3) 

(3.5) 

(3.4) 

(1.9) 

F 

5 

1 

< 0.1 

0.1 

- 

u 

D 

(44) 

3 

1.6 

0.6 

1.5 

E 

(49) 

(8) 

(1.9) 

(1.3) 

(1.5) 

F 

13 

2 

1.3 

0.7 

1.0 
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Table 3 Results of a 15-parameters fit to the HERA I data and to the HERA I and Tevatron W- and Z-boson data. The 
contribution to the total Xmin each data set and the corresponding number of points are shown. 




HERA I 

HERA I -H Tevatron W, Z 

Data set 


X^ / number of points 

X^ / number of points 

NC DIS cross sections Hl-ZEUS combined 

e~p. 

112 / 145 

106 / 145 

NC DIS cross sections Hl-ZEUS combined 

e+p. 

326 / 337 

334 / 337 

CC DIS cross sections Hl-ZEUS combined 

e~p. 

20 / 34 

19 / 34 

CC DIS cross sections Hl-ZEUS combined 

e“fp. 

27 / 34 

33 / 34 

HERA I correlated x^ 


21 

23 

DO da{Z)/dy 


- 

23 / 28 

CDF dcj{Z)ldy 


- 

33 / 28 

DO muon charge asymmetry in W fiu 


- 

12 / 10 

CDF W charge asymmetry in W ^ ei/ 


- 

15 / 13 

DO W charge asymmetry in W —>■ ei/ 


- 

16 / 14 

Total xLn / dof 


505 / 535 

615 / 628 


Table 4 Comparison between the Tevatron lE-boson measurements and the CTlOnlo, MMHT2014, and NNPDF3.0 PDFs. 
The partial of each data set, the total x^, and the total x^ without PDFs uncertainties are shown. 


PDF set 

CTlOnlo 

MMHT2014 

NNPDF3.0 


X^ / number of points 

X^ / number of points 

X^ / number of points 

DO muon charge asymmetry in VF —>■ y,v 

13 / 10 

- 

12 / 10 

CDF W charge asymmetry in W ^ eu 

15 / 13 

- 

17 / 13 

DO W charge asymmetry in VF —J- ei/ 

15 / 14 

10/14 

3/14 

PDF correlated x^ 

4 

3 

6 

Total x^ / dof 

47 / 37 

14 / 14 

39 / 37 

Total x^ / dof without PDFs uncertainties 

368/37 

42/14 

846 / 37 


Table 5 Results of a 15-parameters fit to the HERA I and Tevatron VF-boson lepton asymmetry and Z-boson data. The 
contribution to the total Xmin of each data set and the corresponding number of points are shown. 




HERA I + Tevatron W-lepton, Z 

Data set 


X^ / number of points 

NC DIS cross sections Hl-ZEUS combined 

e~p. 

107 / 145 

NC DIS cross sections Hl-ZEUS combined 

e"*" p. 

334 / 337 

CC DIS cross sections Hl-ZEUS combined 

e~p. 

20 / 34 

CC DIS cross sections Hl-ZEUS combined 

e"*" p. 

32 / 34 

HERA I correlated x^ 


22 

DO dij{Z)/dy 


23/ 28 

CDF da{Z)/dy 


32/ 28 

DO muon charge asymmetry in IF —>■ yv 


13 / 10 

DO electron charge asymmetry in W ^ ei' 


19 / 13 

Total Xmin / dof 


603 / 614 


















12 


HERAFitter developers’ team: S. Camarda^ et al. 



a. 


>. 

"d 

■D 


(0 

(0 

Q 


o 

0 > 


120 

100 

80 

60 

40 

20 

0 

-20 

-40 

-60 

1.2 

1 

0.8 


pp ^ Z-> eV; Vs = 1.96 TeV 

-f 

J L = 2.1 fb’ 



CDF Data 
♦ 8 uncorrelated 
8 total 

- Theory 

- - - Theory + shifts 


15-parameters fit 

I I 1_1_1_1 L 
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(b) 


Figure 1 Theoretical predictions evaluated with the PDFs extracted from a fit to the HERA I and Tevatron data are 
compared to (a) Z-boson differential cross section as a function of rapidity, measured by the DO collaboration and (b) Z-boson 
differential cross section as a function of rapidity, measured by the CDF collaboration. The red continuous lines correspond 
to the theoretical predictions, the red dashed lines are the theoretical predictions shifted by the experimental shift terms 

yellow bands show the total experimental uncertainty, the black vertical bars show the quadratic 
sum of statistical and systematic uncorrelated uncertainties. Note, the theoretical predictions and the theoretical predictions 
shifted by the experimental shift terms are nearly identical, and virtually indistinguishable in these plots. 
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h,l 


(a) 




Figure 2 Theoretical predictions evaluated with the PDFs extracted from a fit to the HERA I and Tevatron data are compared 
to (a) the charge asymmetry of muons as a function of rapidity in W —>■ fj,u decays, measured by the DO collaboration, (b) the 
IT-boson charge asymmetry as a function of rapidity in the W eu decay channel, measured by the CDF collaboration, and 
(c) the VF-boson charge asymmetry as a function of rapidity in the W eu decay channel, measured by the DO collaboration. 
The red continuous lines correspond to the theoretical predictions, the red dashed lines are the theoretical predictions shifted 
by the experimental shift terms of Eq. (2). The yellow bands show the total experimental uncertainty, the black 

vertical bars show the quadratic sum of statistical and systematic uncorrelated uncertainties. 
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Figure 3 PDFs at the starting scale = 1.7 GeV^ as a function of Bjorken-a; for (a) u^, (b) (c) u, and (d) d, determined 

with a fit to the HERA 1 data (blue), and with a fit to the HERA 1 and Tevatron W- and Z-boson data (yellow). 
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(c) 


(d) 


Figure 4 Relative PDF uncertainties at the starting scale = 1.7 GeV^ as a function of Bjorken-a; for (a) (b) (c) 

u, and (d) d, determined with a fit to the HERA I data (blue), and with a fit to the HERA I and Tevatron W- and Z-boson 
data (yellow). 
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Figure 5 (a) d-valence PDF at the scale = 1.7 GeV^ as a function of Bjorken-a; and (b) d-valence relative PDF 

uncertainties, determined with a fit to the HERA 1 data (blue), with a fit to the HERA I and Tevatron VF-boson asymmetry 
and Z-boson data (yellow), and with a fit to the HERA 1 and Tevatron VF-boson lepton asymmetry and Z-boson data (green). 
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Figure 6 d-valence PDF at the scale = 1.7 GeV^ as a function of Bjorken-i; before and after profiling for the (a) CTlOnlo 
and (c) MMHT2014 PDFs and the corresponding relative uncertainties for (b) CTlOnlo and (c) MMHT2014. 
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Figure 7 (a) d-valence PDF at the scale = 1.7 GeV^ as a function of Bjorken-r determined from a fit to the HERA 1 

and Tevatron W- and Z-boson data, and from the CTlOnlo and MMHT2014 PDFs; (b) ratio of d-valence PDFs central values 
and uncertainties with respect to the d-valence PDF determined from a fit the HERA I and Tevatron W- and Z-boson data. 




Figure 8 (a) d-valence PDF at the scale = 1.7 GeV^ as a function of Bjorken-r determined from a fit to the HERA 1 and 

Tevatron W- and Z-boson data, and from the profiled GTlOnlo and MMHT2014 PDFs; (b) ratio of d-valence PDFs central 
values and uncertainties with respect to the d-valence PDF determined from a fit the HERA I and Tevatron W- and Z-boson 
data. 












































